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Abstract 

We study the dynamics of iterates at the transition to chaos in 
the logistic map and find that it is constituted by an infinite family 
of Mori's g-phase transitions. Starting from Feigenbaum's a function 
for the diameters ratio, we determine the atypical weak sensitivity 
to initial conditions £t associated to each g-phase transition and find 
that it obeys the form suggested by the Tsallis statistics. The specific 
values of the variable q at which the g-phase transitions take place 
are identified with the specific values for the Tsallis entropic index q 
in the corresponding £t. We describe too the bifurcation gap induced 
by external noise and show that its properties exhibit the character- 
istic elements of glassy dynamics close to vitrification in supercooled 
liquids, e.g. two-step relaxation, aging and a relationship between 
relaxation time and entropy. 

1 Introduction 

The logistic equation was introduced in 1845 by the Belgian mathematician 
and sociologist Pierre-Frangois Verhulst to model the growth of populations 
limited by finite resources [1]. The discrete time variable version of Verhulst's 
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growth law, the logistic map, has become a foundation stone for the theory 
of nonlinear dynamics. The logistic map is the archetypal example of how 
the use as starting point of simple non-linear discrete maps have often led to 
significant developments in the theory of non-linear dynamical systems [2]. 
The discovery of the universal properties associated to the renowned period- 
doubling and intermittency routes to chaos displayed by the logistic map, 
similar to those of conventional critical phenomena in statistical physics, 
triggered, about three decades ago, an upsurge of activity in the field and 
now both routes, as well as many other remarkable features displayed by the 
logistic map, are well understood. 

Here we might argue that at the present time the logistic map is again 
becoming a prototypical model system. This time for the assessment of the 
validity and understanding of the reasons for applicability of the nonextensive 
generalization of the Boltzmann Gibbs (BG) statistical mechanics [3], [4]. 
This is a formalism assumed to be appropriate for circumstances where the 
system is out of the range of validity of the canonical BG theory. And these 
circumstances are believed in some cases to be a breakdown in the chain of 
increasing randomness from non-ergodicity to completely developed phase- 
space mixing. 

The logistic map contains infinite families of critical attractors at which 
the ergodic and mixing properties breakdown. These are the tangent bi- 
furcations and the accumulation point (s) of the pitchfork bifurcations, the 
so-called onset of chaos [2]. At each of the map critical attractors the Lya- 
punov coefficient Ai vanishes, and the sensitivity to initial conditions £ t for 
large iteration time t ceases to obey exponential behavior, exhibiting instead 
power-law or faster than exponential behavior. The pitchfork bifurcations 
are also critical attractors at which the negative Lyapunov coefficient of pe- 
riodic orbits goes to cero. There are other attractors at which the Lyapunov 
coefficient diverges to minus infinity, where there is faster than exponential 
convergence of orbits. These are the superstable attractors located between 
successive pitchfork bifurcations and their accumulation point is also the 
onset of chaos. 

Here we review briefly specific and rigorous results on the dynamics asso- 
ciated to the critical attractors of the logistic map, or of its generalization to 
non-linearity of order z > 1, f^(x) = 1 — (J, \x\ z , — 1 < x < 1, < // < 2. (The 
phase space variable is x, the control parameter is /i and the conventional 
logistic map corresponds to z — 2). Our results relate to the anomalous sen- 
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sitivity to initial conditions at the onset of chaos, the associated spectrum of 
Tsallis g-Lyapunov coefficients, and the relationship of these with the Mori 
g-phase transitions [5], one of which was originally observed numerically for 
the Feigenbaum attractor [5] [6]. In particular, we identify the Mori singu- 
larities in the Lyapunov spectra with the appearance of special values for 
the Tsallis entropic index q. As the properties of the logistic map are very 
familiar and well understood it is of interest to see how previous knowledge 
fits in with the new perspective. 

Tsallis suggested [7] that for critical attractors £ t (defined as £t(xo) = 
\im.Ax ^o(Axt/ Axo) where Axo is the initial separation of two orbits and 
Ax t that at time t), has the form 

&(xo) = expJA^o) t] = [l-(q- l)X q (x ) t]' 1 ^, (1) 

that yields the customary exponential £ t with Ai when q — > 1. In Eq. (1) q is 
the entropic index and X q is the g-generalized Lyapunov coefficient; exp q (x) = 
[1 — (q — l):r] _1// ( 9_1 ) is the g-exponential function. Tsallis also suggested [7] 
that the Pesin identity K\ = Ai (where the rate of entropy production K 1 
is given by Kit = Si(t) — Si(0), t large, and Si = — Pi hi Pi) would be 
generalized to K q = X q , where the g-generalized rate of entropy production 
K q is defined via K q t = S q (t) — S q (0), t large, and where 

s ff = 5>HY-) = l -^A (2 ) 

i \PiJ 5-1 

is the Tsallis entropy; ln q y = (y 1 ~ q — — q) is the inverse of expJy). 

To check on the above, we have analyzed recently [8]- [11] both the pitch- 
fork and tangent bifurcations and the onset of chaos of the logistic map and 
found that indeed the Tsallis suggestions hold for these critical attractors, 
though with some qualifications. For the case of the tangent bifurcation it 
is important to neglect the feedback mechanism into the neighborhood of 
the tangency to avoid the crossover out of the g-exponential form for £ t . 
As we explain below, for the onset of chaos there is a multiplicity of q- 
indexes, appearing in pairs qj, and Qj = 2 — qj, j = 0, 1, and a spectrum 
of g-Lyapunov coefficients A^'^ and Xq^ for each qj and Qj, respectively. 
(The superindex (k,l) refers to starting and final trajectory positions). The 
dynamics of the attractor confers the g-indexes a decreasing order of im- 
portance. Retaining only the dominant indexes go and Qq yields a quite 
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reasonable description of the dynamics and considering the next few leading 
indexes provides a satisfactorily accurate account. 

Although in brief, we also describe our finding [12] that the dynamics at 
the noise-perturbed edge of chaos in logistic maps is analogous to that ob- 
served in supercooled liquids close to vitrification. The three major features 
of glassy dynamics in structural glass formers, two-step relaxation, aging, 
and a relationship between relaxation time and configurational entropy, are 
displayed by orbits with vanishing Lyapunov coefficient. The known prop- 
erties in control-parameter space of the noise-induced bifurcation gap play a 
central role in determining the characteristics of dynamical relaxation at the 
chaos threshold. 

2 Critical attractors in the logistic map 

For our purposes it is convenient to recall some essentials of logistic map 
properties. The accumulation point of the period doublings and also of the 
chaotic band splittings is the Feigenbaum attractor that marks the threshold 
between periodic and chaotic orbits, at /^(z), with = 1.40115... when 
z = 2 . The locations of period doublings, at fj, — /i n < //qo, and band 
splittings, at /i — fi n > /i^, obey for large n power laws of the form fi n —fioo ~ 
8{z)~ n and /ioo — fi n ~ 8{z)~ n , with 5 = 0.46692... when z = 2, which is one 
of the two Feigenbaum's universal constants. For our use below we recall also 
the sequence of parameter values JZ n employed to define the diameters d n of 
the bifurcation forks that form the period-doubling cascade sequence. At \i — 
~p n the map displays a 'superstable' periodic orbit of length 2 n that contains 
the point x — 0. For large n the distances to x — 0, of the iterate positions 
in such 2 n -cycle that are closest to x = 0, d n = ; (0), have constant 

ratios d n /d n+ i = —a(z); a = 2.50290... when z = 2, which is the second 
of the Feigenbaum's constants. A set of diameters with scaling properties 
similar to those of d n can also be defined for the band splitting sequence [2]. 
Other diameters in the 2 n -supercycles are defined as the distance of the mth 
element x m to its nearest neighbor /A (x m ). That is 

dn, m = f^ Vt ~ 1) (0)-fS: ) (0) i m = 0,1,2,..., (3) 
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with d ni o = d n . Feigenbaum [13] constructed the auxiliary function 



a n (m) 



d 



'n+l,m 
dn,m 



(4) 



to quantify the rate of change of the diameters and showed that it has finite 
(jump) discontinuities at all rationals, as can be seen by considering the 
variable y = m/2 n+1 with n large (and therefore omitting the subindex n). 
One obtains [13] [2] a(0) = — 1/cc, but <x(0 + ) = l/a z , and through the 
antiperiodic property a(y + 1/2) = —a(y), also cr (1/2) = 1/a, but o{l/2 + 
+ ) = — l/a z . Other discontinuities in a(y) appear at y = 1/4, 1/8, etc. As 
these decrease rapidly in most cases it is only necessary to consider the first 
few. 

An important factor of our work is that the sensitivity to initial conditions 
£t(#o) can be evaluated for trajectories within the Feigenbaum attractor via 
consideration of the discontinuities of a n (m). If the initial separation is 
chosen to be a diameter Ax = d n>m and the final time t is chosen to have 
the form t = 2 n — 1, then Ax t = dra,m+2"-i, and £t( x o) can De written [14] as 



For clarity of presentation, we shall only use absolute values of positions so 
that the dynamics of iterates do not carry information on the self-similar 
properties of "left" and "right" symbolic dynamic sequences [2]. This choice 
does not affect results on the sensitivity to initial conditions. 

3 Mori's (/-phase transitions in the logistic 
map 

During the late 1980 's Mori and coworkers developed a comprehensive ther- 
modynamic formalism to characterize drastic changes at bifurcations and at 
other singular phenomena in low dimensional maps [5]. The formalism was 
also adapted to the study of critical chaotic attractors and was illustrated by 
considering the specific case of the onset of chaos in the logistic map [5] [6] 
[15]. For critical attactors the scheme involves the evaluation of fluctuations 




large. 
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of the generalized finite-time Lyapunov coefficient 



f * _1 
X(t,x Q ) = t— E ln 

111 L i=Q 



t > f. 



(6) 



Notice the replacement of the customary t by hit above, as the ordinary 
Lyapunov coefficient 



Ai(xo) 



f * _1 
lim - In 



dfu 



[Xi 



(7) 



vanishes at /i^, t — > oo. 

The probability density distribution for the values of A, P(A, £), is written 
in the form 

P(\,t)=t-* ( - x >P(p,t), f»l, (8) 

where ip(X) is a concave spectrum of the fluctuations of A with minimum 
■0(0) = and is obtained as the Legendre transform of the 'free energy' 
function 0(g), defined as 



= - lim r~7 lnZ ( t '?)' 
mi 

where 5) is the dynamic partition function 



Z(t,q) = J dX P(X,t) t-to-V x . 



(9) 



(10) 



The 'coarse-grained' function of generalized Lyapunov coefficients X(q) is 
given by X(q) = d(f)(q)/dq and the variance v(q) of P(A, t) by v(q) = dX(q)/dq 
[5] [6]. Notice the special weight t _ (' ? ~ 1 ) A in the partition function Z(t,q) as 
this shapes the quantities derived from it. These functions are the dynamic 
counterparts of the Renyi dimensions D q and the spectrum f(a) that char- 
acterize the geometric structure of the attractor. 

A " g-phase" transition is indicated by a section of linear slope m c = 1 — q c 
in the spectrum (free energy) ip(X), a discontinuity at q c in the Lyapunov 
function (order parameter) X(q), and a divergence at q c in the variance (sus- 
ceptibility) v(q). For the onset of chaos at /i o(- 2 = 2) a single g-phase tran- 
sition was numerically determined [5] [6] [15] and found to exhibit a value 
close to m c = — (1 — g c ) ~ — 0.7; arguments were provided for this value to 
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be m c = — (1 — q c ) = — ln2/lno: = —0.7555... Our analysis described be- 
low shows that the older results give a broad picture of the dynamics at the 
Feigenbaum attractor and that actually an infinite family of g-phase transi- 
tions of decreasing weights take place at ji^. 



4 Tsallis dynamics at the edge of chaos 

By taking as initial condition xq = we found that the resulting orbit consists 
of trajectories made of intertwined power laws that asymptotically reproduce 
the entire period-doubling cascade that occurs for /x < /j^ [9] [11]. This 
orbit captures the properties of the so-called 'superstable' orbits at ~p n < fi^, 
n = 1,2, ... [2] (see Fig. 1), and can be used as reference to read all other 
orbits within the attractor. At /i^ the Lyapunov coefficient Ai vanishes and 
in its place there appears a spectrum of g-Lyapunov coefficients A^ fc '^. This 
spectrum was originally studied in Refs. [15] and [6] and our recent interest 
has been to study its properties in more detail to examine their relationship 
with the Tsallis statistics. Recent analytical results about the g-Lyapunov 
coefficients and the g-generalized Pesin identity are given in Refs. [9] and 
[11]. 

Now, consider a first approximation to the function a n (m) for n large, 
applicable to general non-linearity z > 1. This is to assume that half of the 
diameters scale as a (as in the most sparse region of the attractor) while the 
other half scale as «o — °- z ( as m the most crowded region of the attractor). 
This approximation captures the effect on £ t of the most dominant trajecto- 
ries within the attractor. With these two scaling factors cr n (m) becomes the 
periodic step function 

a = a z , < m < 2 n -\ 

a, 2 n ~ l < m < 2-2™-\ 

a n (m) ~ ) -a = -a z , 2 • 2 n ~ l < m < 3 • 2"-\ 1 j 

-a , 3-2"- 1 <m<A-2 n '\... 

Use of Eqn. (11) into Eqn. (5) for the sensitivity £t(m) yields the result 



a 



-(*-i)n m = (2k + l)2 n - 1 



6(m) " 1 a^ n ,' m = (2k + 2)2"-\ (12) 
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Figure 1: Absolute values of positions in logarithmic scales of iterations r for 
a trajectory at ji^ with initial condition x = 0. The numbers correspond 
to iteration times. The power-law decay of the time subsequences can be 
clearly appreciated. 



where k = 0, 1, ... and where the final observation time is of the form t = I2 n — 
1, I = 1,2, .... With the introduction of the total time variable r = m + 1 + 1 
Eq. (12) can be rewritten in terms of the g-exponential functions 



f r (m) 



[1 + (1 - ^A^V] 1 /^) ) m = (2fc + 2)2-1, 
[1 + (1 - g )A^ ) r] 1 /d-Qo) ? m = (2k + l)2 n ~\ 



where 



5o = 1 - 



In 2 



1 - 



In 2 



\(k,l) 
90 



lnao/a (z — 1) In a' 

z — 1) In a 



Qo = 1 + 



(fc + 2Z + l)ln2' 
In 2 In 2 



In a /a 



1 + 



'z — 1) In a' 



(13) 

(14) 
(15) 
(16) 
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A 



(k,l) 
Qo 



2(z - Dlnct 



(17) 

(2A; + 4/ + l)ln2 V ; 

Notice that Q = 2 — q . For z = 2 one obtains Qo ~ 1.7555 and g ~ 0.2445, 
this latter value agrees with that obtained in several earlier studies [7] [16] 
[9] [11]. The two scales considered describe correctly only trajectories that 
start at the most sparse region of the multifractal (x ~ 0) and terminate at 
its most crowded region (x ~ 1), or the inverse. (This is why we obtain the 
two conjugate values q and Q = 2 — q, as the inverse of the g-exponential 
function satisfies exp q (y) = 1/ exp 2 _ ? (— y)). The vertical lines in Fig. 2a 

represent the ranges of values obtained when z = 2 for A^'^ and Xq^ ■ See 
Ref. [14] for more details. 
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Figure 2: g-phase transitions with index values go an d Qo — 2 — go obtained from 
the main discontinuity in a n (m). See text for details. 



We consider the next discontinuities of importance in a n (m). Indepen- 
dently of the number of discontinuities taken into account one obtains g- 
exponential forms for £ t . The value of a(y — m/2 n+1 ) at y — 1/4 measures 
1/cti, with cti ~ 5.4588 for z = 2, and it is associated to one 'midway' 
region between the most crowded and most sparse regions of the attractor 
(the other 'midway' region being associated to a (3/4)). With three scaling 
factors, a, ao and a±, we have three values for the q index, go, gi and g2 
(together with the conjugate values Qo = 2 — g , Q\ = 2 — gi and Q2 = 2 — g 2 
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for the inverse trajectories). To each value of g there is a set of g-Lyapunov 
coefficients running from a maximum A giimax to zero (or a minimum AQ iim i n 
to zero). The results when z = 2 for the ranges of values obtained for the 
g-Lyapunov coefficients are shown as the vertical lines in Figs. 2a, 3a and 
4a. Similar results are obtained for the case of four discontinuities in a n (m), 
etc. 

5 A family of (/-phase transitions at the edge 
of chaos 

As a function of the variable — oo < q < oo the g-Lyapunov coefficients ob- 
tained in the previous section are functions with two steps with the jumps 
located at q = qi = 1 — ln2/ In oti(z) / a(z) and q = Qt = 2 — q^. Imme- 
diate contact can be established with the formalism developed by Mori and 
coworkers and the q phase transition obtained in Ref . [6] . Each step function 
for \{q) can be integrated to obtain the spectrum 0(g) (X(q) = d(f>/dq) and 
from this its Legendre transform (= — (1 — g)A). We illustrate this 
first with a n (m) obtained with two scale factors, as in Eq. (11). We show 
numerical values for the case z — 2. From Eqs. (14) to (17) we obtain 



where \ qo , max = In a/ In 2 ~ 1.323 and AQ 0jmin = —2 In a/ In 2 ~ —2.646. 

The free energy functions 0(g) and ^(A) that correspond to Eq. (18) are 
given by 



We show these functions in Fig. 2. The constant slopes of represent 
the g-phase transitions associated to trajectories linking two regions of the 
attractor, that in this case are its most crowded and most sparse, and their 



A(?) 



' A 30imax , q < q = 1 - In 2/(z - 1) In a ~ 0.2445, 
< 0, qo<q< Qo, 

K AQ ,min, q > Qo = 2 - g - 1-7555, 



(18) 




and 
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values 1 — go and go — 1 correspond to those obtained for the g-exponential £t 
Eq. (13). The slope g — 1 ~ —0.7555 coincides with that originally detected 
by Mori and colleagues [5]. 

When we consider also the next discontinuity of importance in a n (m), at 
<j (1/4) = 1/cci, we obtain a set of two g-phase transitions for each of the 
three values of the q index, g , q 1 and q^. We show in Figs. 2, 3 and 4 the 
functions X(q), 0(g), and ip(\) obtained for these three cases. The parameter 
values for the g-phase transitions at 1 — go and go — 1 appear again, but now 
we have also two other sets at 1 — qi and q± — 1, and at 1 — g2 and g2 — 1, 
that correspond, respectively, to orbits that link a 'midway' region of the 
attractor with the most sparse region, and with the most crowded region of 
the attractor. 
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Figure 3: g-phase transitions with index values gi and Q\ = 1 — q\ obtained from 
the main two discontinuities in a n (m). See text for details. 



6 Noisy dynamics at the edge of chaos 

Consider now the logistic map z = 2 in the presence of additive noise 

x t+ i = U(x t ) = l-nx 2 t + Xt e, -1 < x t < 1, < n < 2, (19) 
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Figure 4: g-phase transitions with index values q2 and Q2 = 2 — q2 obtained from 
the main two discontinuities in a n (m). See text for details. 



where Xt is the random variable with average (xtXf) — $t.t', and e measures 
the noise intensity [2] [17]. Except for a set of zero measure, all the tra- 
jectories with fioo(e = 0) and initial condition — 1 < x$ < 1 fall into the 
attractor with fractal dimension df = 0.5338.... These trajectories represent 
nonergodic states, since as t — > 00 only a Cantor set of positions is acces- 
sible out of the total phase space. For e > the noise fluctuations erase 
the fine features of the periodic attractors as these widen into bands similar 
to those in the chaotic attractors, yet there remains a definite transition to 
chaos at /ioo(^) where the Lyapunov exponent Ai changes sign. The period 
doubling of bands ends at a finite value as the edge of chaos transition 
is approached and then decreases at the other side of the transition. This 
effect displays scaling features and is referred to as the bifurcation gap [2] 
[17]. When e > the trajectories visit sequentially a set of 2™ disjoint bands 
or segments leading to a cycle, but the behavior inside each band is com- 
pletely chaotic. These trajectories represent ergodic states as the accessible 
positions have a fractal dimension equal to the dimension of phase space. 
Thus the removal of the noise e — > at leads to an ergodic to nonergodic 
transition in the map. 

In the absence of noise (e = 0) the diameter positions £2™ = d n = a~ n 
visited at times r = 2 n by the trajectory starting at xq — 1 is given by the 
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Feigenbaum fixed-point map solution g(x), 



9 (&o) 



--1/1-90 



(20) 



that in turn is obtained from the n — > oo convergence of the 2 n th map 
composition to (— a)~ n g(a: n £) with ct = 2 1 /( 1 ~' ?0 ). When x = one obtains 
in general [9] 



/ (2J+l) (0 ^(2»-i) (0) 



,(21+1), 



(0) a~ n ,T= (2Z + l)2 n , /,n = 0,1,... 



(21) 

When the noise is turned on (e always small) the 2 n th map composition 
converges instead to 



'-a)- n [g(a n x) + X eK n G A (a n x)}, 



(22) 



where k a constant whose numerically determined [18], [19] value k ~ 6.619 
is well approximated by v = 2-\/2a(l + l/a 2 )^ 1 / 2 , the ratio of the intensity of 
successive subharmonics in the map power spectrum [19], [2]. The connection 
between k and the e-independent v stems from the necessary coincidence of 
two ratios, that of noise levels causing band-merging transitions for succes- 
sive 2™ and 2 n+1 periods and that of spectral peaks at the corresponding 
parameter values ji n and /i n +i [19], [2]. Following the same procedure as 
above we see that the orbits x T at fiooi^) satisfy, in place of Eq. (20), the 
relation 



X T = T 



-1/1-90 



gir^-^x) + X er 1/1 - r G A (T 1/1 - qo x) 



(23) 



where G\(x) is the first-order perturbation eigenfunction, and where r 
1 — In 2/ In k ~ 0.6332. So that use of x = yields 



-1/1-50 



1+X£T 



1/1-r 



(24) 



or 

x t = exp 2 _ qo (-\ qo t) [1 + x^exp r (A r t)] (25) 

where t — r — 1, X qo = In a/ In 2 (A ?0;max of the previous section) and A r = 
ln/t/ln2. 

At each noise level e there is a 'crossover' or 'relaxation' time t x = r x — 1 
when the fluctuations start destroying the detailed structure imprinted by 
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the attractor on the orbits with xq = 0. This time is given by r x = e r_1 , 
the time when the fluctuation term in the perturbation expression for x T 
becomes e-independent, i.e. 

Zr^T-^ll + xl- (26) 

Thus, there are two regimes for time evolution at //00(e). When t < t x the 
fluctuations are smaller than the distances between adjacent subsequence 
positions of the noiseless orbit at /^(O), and the iterate positions in the 
presence of noise fall within small non overlapping bands each around the 
e = position for that r. In this regime the dynamics follows in effect the 
same subsequence pattern as in the noiseless case. When r ~ t x the width 
of the fluctuation-generated band visited at time r x = 2 N matches the dis- 
tance between two consecutive diameters, — d^+i where N ~ — In ej In k, 
and this signals a cutoff in the advance through the position subsequences. 
At longer times r > r x the orbits are unable to stick to the fine period- 
doubling structure of the attractor. In this 2nd regime the iterate follows 
an increasingly chaotic trajectory as bands merge progressively. This is the 
dynamical image - observed along the time evolution for the orbits of a single 
state /too( £ ) - °f the static bifurcation gap first described in the map space 
of position x and control parameter /t [17], [18]. 



7 Analogy with glassy dynamics 

We recall the main dynamical properties displayed by supercooled liquids on 
approach to glass formation. One is the growth of a plateau and for that 
reason a two-step process of relaxation, as presented by the time evolution of 
correlations e.g. the intermediate scattering function F k [20]. This consists 
of a primary power-law decay in time t (so-called relaxation) that leads 
into the plateau, the duration t x of which diverges also as a power law of 
the difference T — T g as the temperature T decreases to a critical value 
T g . After t x there is a secondary power law decay (so-called 'a' relaxation) 
that leads to a conventional equilibrium state [20]. A second characteristic 
dynamical property of glasses is the loss of time translation invariance, a 
characteristic known as aging [21]. The time decay of relaxation functions 
and correlations display a scaling dependence on the ratio t/t w where t w 
is a waiting time. A third distinctive property is that the experimentally 
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observed relaxation behavior of supercooled liquids is described, via regular 
heat capacity assumptions [20], by the so-called Adam-Gibbs equation, 



t x = Aexp(B/TS c ) 



(27) 



where the relaxation time t x can be identified with the viscosity, and the 
configurational entropy S c is related to the number of minima of the fluid's 
potential energy surface (and A and B are constants). 

Returning to the map, phase space is sampled at noise level e by orbits 
that visit points within the set of 2 N bands of widths ~ e, and this takes place 
in time in the same way that period doubling and band merging proceeds in 
the presence of a bifurcation gap when the control parameter is run through 
the interval < /i < 2. That is, the trajectories starting at x = duplicate 
the number of visited bands at times r = 2 n , n = 1, ...,N, the bifurcation 
gap is reached at r x = 2 N , after which the orbits fall within bands that 
merge by pairs at times r = 2 N+n , n = 1,...,N. The sensitivity to initial 
conditions grows as £ t = exp (X qo t) (q — 1 — ln2/lnct < 1) for t < t x , 
but for t > t x the fluctuations dominate and £ t grows exponentially as the 
trajectory has become chaotic and so one anticipates an exponential £ f (or 
q — 1). We interpret this behavior to be the dynamical system analog of 
the 'a' relaxation in supercooled fluids. The plateau duration t x — ■> oo as 
e — > 0. Additionally, trajectories with initial conditions Xq not belonging to 
the attractor exhibit an initial relaxation stretch towards the plateau as the 
orbit falls into the attractor. This appears as the analog of the relaxation 
in supercooled liquids. See [12]. 

Next, we determine the entropy of the orbits starting at xq = as they 
enter the bifurcation gap at t x (e) when the maximum number 2^ of bands 
allowed by the fluctuations is reached. The entropy S c (/j,oo, 2 N ) associated to 
the state at //«>(£) has the form 



and where p(x) is the distribution for the noise random variable. In terms 
of t x , given that 2^ = 1 + t x and e — (1 + t x )~ 1 / 1 ~ r , one has 



5^00,2^ = 2^, 
since each of the 2^ bands contributes with an entropy es where 



(28) 




(29) 



5 C (^oo,^)/s=(l+t,)- r/1 - r 



(30) 
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or, conversely, 

t x = (s/S c )^/ r . (31) 

Since t x ~ £ r_1 , r — 1 ~ —0.3668 and (1 — r)/r ~ 0.5792 then ^ — > oo and 
— > as e — > 0, i.e. the relaxation time diverges as the 'landscape' entropy- 
vanishes. We interpret this relationship between t x and the entropy S c to be 
the dynamical system analog of the Adam-Gibbs formula for a supercooled 
liquid. See [12]. Notice that Eq.(31) is a power law in S*" 1 while for structural 
glasses it is an exponential in S' 1 [20]. 

Last, we examine the aging scaling property of the trajectories x t at 
fioo(e). The case e = is more readily appraised because this property 
is actually built into the position subsequences x T = g^ T \0) , r = (21 + 
l)2 n , l,n — 0,1,... These subsequences are relevant for the description of 
trajectories that are 'detained' at a given attractor position for a waiting 
period of time t w and then 'released' to the normal iterative procedure. We 
chose the holding positions to be any of those along the top band shown in 
Fig. 1 for a waiting time t w = 21 + 1, I = 0, 1, .... Notice that, as shown 
in Fig. 1, for the x = orbit these positions are visited at odd iteration 
times. The lower-bound positions for these trajectories are given by those of 
the subsequences at times {21 +1)2™ (see Fig. 1). Writing r as r = t w + 1 we 
have that t/t w = 2 n — 1 and 

xt +t ,=9 M (0)9 {tM (0) (32) 

or 

x t+tw = g^\0) ex P J-\ qo t/t w ). (33) 

This property is gradually modified when noise is turned on. The presence 
of a bifurcation gap limits its range of validity to total times t w + 1 < t x (e) 
and so progressively disappears as e is increased. See Ref. [12]. 



8 Concluding remarks 

We have re-examined the dynamical properties at the onset of chaos in the 
logistic map and obtained further understanding about their nature. We 
exhibited links between original developments, such as Feigenbaum's a func- 
tion, Mori's g-phase transitions and the noise-induced bifurcation gap, with 
more recent advances, such as g-exponential sensitivity to initial conditions 
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[9], g-generalized Pesin identity [11], and dynamics of glass formation [12]. 
The dynamics at the edge of chaos is anomalous because it is an incipient 
chaotic attractor with vanishing ordinary Lyapunov coefficient Ai. Chaotic 
orbits possess a time irreversible property that stems from mixing in phase 
space and loss of memory, but orbits within critical attractors are non-mixing 
and have no loss of memory. As a classic illustration of the latter case the 
attractor at the onset of chaos presents dynamical properties with self-similar 
structure that result in a set of power laws for the sensitivity to initial con- 
ditions. We determined exact analytical expressions for £ t . 

Our most striking finding is that the dynamics at the onset of chaos is 
constituted by an infinite family of Mori's g-phase transitions, each associated 
to orbits that have common starting and finishing positions located at specific 
regions of the attractor. Each of these transitions is related to a discontinuity 
in the a function of 'diameter ratios', and this in turn implies a g-exponential 
£ t and a spectrum of g-Lyapunov coefficients for each set of orbits. The 
transitions come in pairs with specific conjugate indexes q and Q = 2 — q, as 
these correspond to switching starting and finishing orbital positions. Since 
the amplitude of the discontinuities in a diminishes rapidly, in practical terms 
there is only need of evaluation for the first few of them. The dominant 
discontinuity is associated to the most crowded and sparse regions of the 
attractor and this alone provides a very reasonable description, as found in 
earlier studies [7] [16] [9] [11]. Thus, the special values for the Tsallis entropic 
index q in £ t are equal to the special values of the variable q in the formalism 
of Mori and colleagues at which the g-phase transitions take place. 

As described, the dynamics of noise-perturbed logistic maps at the chaos 
threshold presents the characteristic features of glassy dynamics observed 
in supercooled liquids. In particular our results are [12]: i) The two-step 
relaxation that takes place when e — > is obtained in terms of the bifurcation 
gap properties, specifically, the plateau duration t x is given by a power law 
in the noise amplitude e. ii) The map analog of the Adam-Gibbs law is given 
also as a power-law relation between t x (e) and the entropy S c (e) associated 
to the noise widening of chaotic bands, iii) The trajectories at //oo(s — * 0) are 
shown to obey a scaling property, characteristic of aging in glassy dynamics, 
of the form x t +t w = h(t w )h( t/t w ) where t w is a waiting time. 

The limit of vanishing noise amplitude e — > (the counterpart of the 
limit T — T g — > in the supercooled liquid) brings about loss of ergodicity. 
This nonergodic state with Ai = corresponds to the limiting state, e — > 
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0, t x — > oo, for a family of small e states with glassy properties, which 
are expressed for t < t x via the g-exponentials of the Tsallis formalism. It 
has been suggested on several occasions [22] [23] that the setting in which 
nonextensive statistics appears to come out is linked to the prevalence of 
nonuniform convergence, such as that involving the thermodynamic iV — > oo 
and the infinitely large time t — > oo limits. Here a similar situation happens, 
that is, if e — > is taken before t — > oo a nonergodic orbit restrained to the 
Feigenbaum attractor and with fully-developed glassy properties is obtained, 
whereas if t — > oo is taken before £->0a chaotic orbit with q — 1 would be 
observed. 
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